Jamming in granular hopper flow 
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Large-scale three dimensional molecular dynamics simulations of hopper flow are presented. The 
flow rate of the system is controlled by the width of the aperture at the bottom. As the steady-state 
flow rate is reduced, the force distribution P(f) changes only slightly, while there is a large change 
in the impulse distribution P(i). In both cases, the distributions show an increase in small forces 
or impulses as the systems approach jamming, the opposite of that seen in previous Lennard- Jones 
simulations. This occurs dynamically as well for a hopper that transitions from a flowing to a 
jammed state over time. The final jammed P(f) is quite distinct from a poured packing P(f) in the 
same geometry. The change in P(i) is a much stronger indicator of the approach to jamming. The 
formation of a peak or plateau in P(f) at the average force is not a general feature of the approach 
to jamming. 



I. INTRODUCTION 

The vertical pressure of grains in a silo becomes inde- 
pendent of depth for sufficiently tall silos [l|, |2|, |3j] . The 
weight of sand in the column is ultimately supported by 
the lateral containing walls. Experiments 0, 0, IE 0' 
theory @, and simulations d floL ITTI | all show that 
these stresses are transmitted in an inhomogeneous man- 
ner. One measure of the inhomogeneity in a granu- 
lar pile is expressed through the distribution of nor- 
mal forces, P(f), where / = F/(F) is the normalized 
force, F is the normal force, and (F) is the average 
force. P(f) of static granular piles exhibits an expo- 
nential tail at lar ge f orces and a plateau or turnover at 
small forces H 1 M El El El 111 E3 • A fundamen- 
tal question is whether these inhomogeneities are present 
in flowing granular material. In addition, if these in- 
homogeneities are distinct in the flowing state, then it 
may be possible to predict the approach of the quiescent 
or "jammed" state by analyzing the change in inhomo- 
geneities as one moves from one to the other. In the 
context of a general conception of jamming |l7j . it has 
been conjectured [18j that the P(f) of a flowing system 
exhibits a decrease at small forces (/ < 1) as the system 
moves from a flowing to a jammed state. 

An ideal system to probe this transition is the hop- 
per E3- Its jammed state (a silo) has been extensively 
studied. For a granular packing of sufficient height, the 
pressure in the depth of the packing is independent of the 
height of the granular matter above it p| . A granular silo 
is also constrained and well-controlled — there is no free 
surface aside from the top. A silo transitions to a hopper 
when an aperture at the base is opened. Numerous exper- 
iments and simulations have examined the transition to 
jamming with the goal of predicting its onset pfl l2ll |22| . 



Here we focus specifically on the distribution of forces and 
impulses in the approach to jamming. 



In experiments the forces between particles are diffi- 
cult to measure directly |23j, so Longhi et al |2J] mea- 
sured the distribution of impulses P(i) (where the nor- 
malized impulse is i = 7/(7), 7 is the impulse, and (7) 
is the average impulse) between particles and the wall 
and related this distribution to the macroscopic behav- 
ior of the system as the system approached jamming. 
This quasi-27) experiment found that the impulse has an 
exponential tail at all flow rates and that at small flow 
rates, the distribution of P(i) increases with decreasing 
flow rate. The distribution of times between collisions 
tends to a power law P(t c ) rj t c 3 ^ 2 , that is, the mean 
time interval tends to diverge just as in the glass tran- 
sition |24j|. Ferguson et al found using 2D event-driven 
frictionless simulations pEj that P(i) exhibits an increase 
in small impulses (i < 1) and argue that this arises from 
the contribution of "frequently-colliding" particles that 
are spatially correlated into 17) collapse strings. 
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Here we present fully 37) molecular dynamics simula- 
tions of hopper flow. We measure the distributions of 
normal forces P{f) and impulses P(i) for the system as 
the system approaches jamming and relate our results to 
previous experiments and 2D simulations. We find that 
P(f) is not a strong indicator of jamming, and that it 
exhibits an increase at small forces as the system ap- 
proaches jamming, exactly the opposite of the predicted 
behavior from earlier simulations 18]. P(f) of a jammed 
system is markedly different from a poured packing in a 
similar geometry, showing that P(f) can depend strongly 
on the history of the system. We also find that P(i) ex- 
hibits an increase at small impulses, in a gree ment with 
previous experiments [24| and simulations |25| . 
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II. SIMULATION METHOD 

Our geometry was inspired by that used in the quasi- 
2D experiments: a conical hopper [2(j. The system is 
periodic in the z direction: particles that fall out of 
the bottom of the system reappear at the top and re- 
fill the container, creating a steady-state flowing system. 
The system is shown in Figure ^a). The top of the 
hopper is a cylindrical container with radius R = lOd, 
where d is the particle diameter. At the bottom of this 
cylindrical region, a cone is joined to the cylinder, with 
radius varying smoothly from R to rj over a distance 
Zf = 50c?. The angle 9 between vertical and the fun- 
nel is 8 = tan -1 f~~^p*Y This system tends to jam for 
rf < 1.5(1 We deliberately make the cylindrical region 
very deep to assure that in the static case the pressure 
is independent of depth. Those particles that leave the 
hopper opening and then rain down on the top of the pile 
quickly lose their kinetic energy and do not affect the flow 
through the opening. The system is prepared by pour- 
ing particles into the top of the container with aperture 
radius r/ = d and a plate at the bottom to prevent the 
exit of particles. After the packing is formed, the plate is 
removed and the aperture widened to the desired radius. 
All measurements are taken after the system has cycled 
particles through at least once and is in a steady state 
unless specifically noted otherwise. We focus our atten- 
tion on the funnel region, and unless otherwise stated, all 
measurements are only of particles in that region. 

We use a molecular dynamics code with a Hertzian 
force law that is described in detail elsewhere 0, |27j • All 
results will be presented in dimensionless units based on 
the mass of a particle m, the diameter of a particle d, 
and the force of gravity g. In this case, the simulations 
used N = 40000 monodis perse particles of diameter d. 
The unit of time is r = y/d/g, the time it takes a par- 
ticle to fall its radius from rest under gravity. Most of 
the simulations were carried out with a spring constant 
k n = 2 x I0 5 mg/d, for which the simulation time step 
St = 10~ 4 t. The particle-particle friction and particle- 
wall friction are the same: /i = fi w = 0.5. 

The flow profile is dependent on the geometry of the 
system, as shown in Fig. [2J The flow is essentially plug- 
like and constant in the cylindrical region of the hopper 
for particles more than one particle away from the wall. 
A one-particle thick layer near the wall flows more slowly 
due to friction with the wall. This plug breaks up in the 
funnel region, with particles in the center of the flow 
flowing more quickly and accelerating as they approach 
the hopper opening. 

Experimental measurements have shown that hoppers 
produce a mass flow rate dependence of M = (r/ — 
wd) 5 / 2 , where w is a constant related to the shape of the 
particles |28| . We see a similar dependence in our simu- 
lations, with w = 1.2, consistent with the experimental 
measurement of w « 1.5 1281. 



III. FORCES AND IMPULSES 

The distribution of particle-particle normal forces P(f) 
for the funnel region is shown in Fig. [3] for a variety of 
aperture radii r/. P(f) is exponential at larg e forces, 
just as observed in static piles 0, ||| 0, 0, hj. All the 
flowing states of the system have similar distributions. 
In Fig. (inset) we show a closeup of P(f) for a range 
of flowing states. There is a slight change in P(f) as 
the system approaches jamming — it displays a slight 
increase at small forces. The peak of the distribution 
moves to smaller forces as the system approaches jam- 
ming. This behavior is the opposite of the behavior seen 
in 2D Lennard- Jones simulations ^^|. This change in 
P(f) is robust and suggests that the behavior of P(f) in 
a jammed system is dependent on the specific physics of 
the interactions between particles. 

We also analyzed the dynamical nature of P(f) for 
rf = lAd. In this case, the system flows freely for a 
short period and then jams. In Fig.0]P(/) is shown for 
this system as a function of time. As the jam forms, the 
distribution abruptly increases at small forces 29] . This 
behavior is compared to that of a poured system with 
rj = d that never flows after the aperture is opened. In 
that case, P(f) is lower at small forces than the flowing 
state. The formation of a jam is history-dependent, and 
the static packing that results from a jam is fundamen- 
tally distinct from a poured packing. P(f) in a jammed 
packing is strongly dependent on the dynamic processes 
used to create the jam and is not equivalent to the P(f) 
of a granular packing created with an alternate method. 
There is no general P(f) that one should expect for all 
jammed packings. 

This hysteretic sudden change in P(f) has also been 
seen in reverse in shearing experiments in cylinders. In 
that case, P(f) is essentially unchanged until yield stress 
is attained, at which point P(f) changes discontinuously 
to a new form [30|. 

We can also probe the dependence of P(f) on the 
particle interactions. We carried out additional simu- 
lations with a much softer normal spring constant, k n — 
2 x 10 3 mg/d. Experiments have shown that very soft par- 
ticle interactions can have a strong effect on P(f) [l2| . 
We compare the behavior of P(f) for the two systems in 
Figure|SJ The trend of increasing the P(f) at small forces 
is the same in both systems, but the increase is more pro- 
nounced and more localized at small forces for the stiffer 
springs. P(f) is thus sensitive to the compressibility of 
the particles and softness of the particle interactions. 

These analyses provides a clue to the discrepancy be- 
tween our simulations and those of O'Hern et al [18( and 
others. Two possible effects separate our simulations and 
theirs. The first and more important effect arises from 
the ahistorical nature of the O'Hern simulations, in which 
the packings are produced by conjugate gradient meth- 
ods, which have no information on the dynamics of the 
transition from flowing state to jammed state. As demon- 
strated above, in granular materials the history of the 
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FIG. 1: (color online) (a) Geometry of the system. The total height of the system is 160d, with the freefall region roughly 20d 
in height, the cylindrical region 90d, and the funnel region 50d. The bottom aperture radius 77 is varied to control the flow. In 
this case, r/ = 3d. All particles with z positions 120d < z < 130d are colored gray (green), (b) Only the tracer particles after 
time t = 300t. (c) Tracer particles only at t = 580r. (d) Tracer particles only at t = 680r. 
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FIG. 2: Velocity profiles for v z at a number of heights z in 
a flowing hopper with r/ = 3d. The profile with the largest 
magnitude was taken at z = Wd, and at every 5d afterward, 
up to a maximum of 130d. \v z | decreases as the funnel widens 
with increasing height, until one reaches the cylindrical re- 
gion, where it is roughly constant. 



system cannot be ignored, and the P(f) of a jammed 
state is quite distinct from that of a poured state. The 
P(f) obtained for the repulsive interactions by O'Hern et 
al more closely resemble those for static granular pack- 
ings and are in fact produced by the quasi-static method 
of conjugate gradient. In our simulations a poured pack- 




f=F/<F> 

FIG. 3: Distribution of normal forces P(f) for flowing hopper 
systems with various apertures r/. P{f) increases slightly 
at small forces as the system approaches jamming. Inset: 
closeup of P(f) at small forces. The lines are guides to the 
eye. 



ing does exhibit a deficit at small forces in relation to 
flowing systems, but all "jammed" systems clearly do 
not exhibit the same characteristics. A poured packing 
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FIG. 4: Distribution of normal forces P(f) for flowing hop- 
per systems with aperture radius 77 = 1.4ci as a function of 
time. The system flow freely until t ~ 3.6r, when a jam devel- 
ops and shifts the system towards an increase in P(f). The 
jammed P(f) is distinct from that of a poured state P(f), 
here shown for a system that never flows after the aperture is 
opened at rj — d. Inset: closeup of P(f) at small forces. 



FIG. 5: Distribution of normal forces P(f) for flowing hopper 
systems with various aperture radius r/. As in the previous 
plots, the symbols represent aperture radii: O r / = l-5d, 
□ 77 = 2d, O 77 = 3d, and < r/ = 5d. Filled-in symbols 
represent stiff springs of k n = 2 x 10 5 mg/d and open symbols 
represent loose springs of k n = 2 x 10 3 mg/ci. Soft springs 
show a smaller effect in the approach to jamming. 



is not jammed in the same sense as a plugged hopper flow 
packing. 

The second effect is that some simulations use an ex- 
tremely soft potential that allows significant overlap. As 
shown above, softer potentials tend to diminish the in- 
crease at small forces observed in our simulations. For 
sufficiently soft potentials, P(f) could have a deficit at 
small forces. Thus, we do not expect our findings to 
agree with simulations undertaken with extremely soft 
potentials. 

To make contact with experiment |24| and simula- 
tion [25l I^J, we also analyze the impulse distribution. 
Figure shows the distribution of particle-particle nor- 
mal impulses P(i) for contacts that terminate in the hop- 
per region . The tails of the system are no l onge r ex- 
ponential; they decline faster than exponential [33] [34| . 
There is also a strong signal as the system approaches 
the jammed state: P{i) increases at small impulses 
(i = 1/(1) < 1). This increase at small impulses is also 
seen in the quasi-2D hopper experiments [2J, [34J . 

Our measurements were begun after the particles had 
already circulated through the container at least once, 
and impulses were calculated for a duration of 25r, which 
is greater in duration than essentially all contacts. Im- 
pulses were calculated for all those contacts that ended in 
the funnel region. These results are somewhat consistent 
with the change of P(i) from 2D event-driven simula- 
tions [25|. In the 2D simulation, the change in P(i) at 
small impulses was controlled by "rapidly colliding par- 
ticles" , which formed ID linear "collapse strings" . We 
denote "rapidly colliding particles" in our case as those 
with the largest number of collisions over a set time win- 
dow, but we observe no obvious spatial correlation of 




i = I/<I> 

FIG. 6: Probability distribution of impulses P(i„) for the 
funnel region. The system exhibits a strong increase in small 
impulses as the system approaches jamming. Inset: closeup 
of P(in) at small impulses. 

these particles. 



IV. CONCLUSIONS 

In summary, our simulations show that for a granular 
hopper system with history-dependent interactions and 
friction, P(f) exhibits an increase in small forces (/ < 1) 
as the system approaches jamming. This increase is ro- 
bust and hysteretic; the P(f) of a jammed packing is dis- 
tinct from the P(f) of a poured packing. The behavior 
of P(f) is somewhat dependent on the interaction: softer 



5 



particle interactions tend to diminish this effect. The his- 
tory dependence of granular materials distin guis hes our 
simulations from many 2D simulations et al |lSj . which 
have no memory of dynamics, and thus predict a gen- 
eral behavior of P(f) that is exactly the opposite of that 
observed in our simulations. The trend is the same for 
the distribution of impulses P(i), which see an increase 
in small impulses (i < 1) as the system approaches jam- 
ming. This trend is also observed in experiments and 
quasi-2D simulations. The behavior of P(f) does not 
appear to be a general feature of jamming, but instead 
depends on the particle interactions and the hysteretic 
quality of the system. 
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